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We argue that all homogeneous and isotropic solutions in nonlinear massive gravity are unstable. 
For this purpose, we study the propagating modes on a Bianchi type-I manifold. We analyze their 
kinetic terms and dispersion relations as the background manifold approaches the homogeneous and 
isotropic limit. We show that in this limit, at least one ghost always exists and that its frequency 
tends to vanish for large scales, meaning that it cannot be integrated out from the low energy effective 
theory. This ghost mode is interpreted as a leading nonlinear perturbation around a homogeneous 
and isotropic background. 



Introduction. The concept of the mass has been 
central in many areas of physics. Gravitation is not an 
exception, and it is one of the simplest but yet unan- 
swered questions whether the graviton, a spin-2 particle 
that mediates gravity, can have a nonvanishing mass or 
not. This question is relevant not only from a theoret- 
ical but also from a phenomenological viewpoint, since 
a nonzero graviton mass may lead to late-time acceler- 
ation of the universe and thus may be considered as an 
alternative to dark energy. 

Recently Refs.0, proposed the first example of a 
fully nonlinear massive gravity theory, where the so called 
Boulware-Deser (BD) ghost 0], which had been one of 
the major obstacles against a stable nonlinear gravity 
theory with a nonvanishing graviton mass, is removed by 
construction. Due to the theoretical and phenomenolog- 
ical motivations mentioned above, this theory has been 
attracting significant interest. 

The first homogeneous and isotropic Friedmann- 
Lemaitre-Robertson- Walker (FLRW) solution for this 
theory was presented in Q for Minkowski fiducial metric 
and then extended to more general fiducial metrics in [f| . 
The analysis of linear perturbations in this general setup 
was carried out also in Although a massive spin-2 
particle generically has 5 propagating degrees of freedom, 
it was found that the number of degrees of freedom in the 
gravity sector was 2, same as in general relativity (GR). 
This is due to the vanishing of the kinetic terms for the 
expected additional degrees 0- This feature may extend 
to other setups: Ref. @ obtained a vanishing kinetic term 
on spherically symmetric inhomogeneous backgrounds. 

The goal of the present paper is to determine the fate 



of the extra degrees of freedom. We find that in the 
nonlinear massive gravity, all cosmological solutions that 
respect homogeneity and isotropy have a ghost, i.e. an 
excitation with a wrong sign kinetic term. Therefore, 
the universe in this theory should have either inhomo- 
geneities [1] or anisotropy [§]. We note that the ghost 
mode found in the present paper is among five degrees of 
freedom of the massive spin-2 field and thus, is not the 
BD ghost. 

The model and the background. Imposing the 
absence of the BD ghost, the massive gravity action, in 
vacuum, can be constructed as 
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1 Backgrounds with additional symmetries which remove the extra 
degrees were introduced in 0. 



with £mg = £2 + 0:3 £3 + 0:4 £4, where 

£ 2 = ±([/C] 2 -[/C 2 ]), 

£ 3 = ±([/C] 3 -3[/C][/C 2 ] + 2[/C 3 ]), 

£4 = 2i([/C] 4 - 6[/C] 2 [/C 2 ] + 3[/C 2 ] 2 + 8[/C][/C 3 ] - 6[/C 4 ]), 

the square brackets denote the trace operation, and 

W^^v-Wo^fYv. (2) 

Here, g^ v and f^ v are physical and fiducial metrics, re- 
spectively. Since we are interested in the stability of the 
gravity sector only, it is sufficient to consider a vacuum 
configuration, with a cosmological constant A. 

The physical metric is chosen to be the simplest 
anisotropic extension of FLRW, namely, the axisymmet- 
ric Bianchi type-I metric 

ds 2 = -N 2 dt 2 + a 2 (e i ' J dx 2 + e~ 2<T Sy dy l dy J ) , (3) 

where N, a, and a are functions of the time variable t. In 
the rest of the paper, Greek indices span the space-time 
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coordinates, while the indices i, j = 2, 3 correspond to the 



coordinates on the y-z plane, with y 



y, V 



z. Since 



our goal is to obtain the stability conditions of this metric 
in the isotropic limit, the whole system in this limit needs 
to reduce to the general cosmological solutions given in 
d, [B| . For this reason, we consider a fiducial metric to 
be in the flat FLRW form, 

ffiv = -rfd^d^ + a^d^d^+Sijd^d^), (4) 



where both n and a are functions of the timc-Stiickelberg 
field cb°. 



The equations of motion for the background can be 
calculated by varying the action with respect to the 
Stiickelberg fields and the metric. As a result, we ob- 
tain three independent equations as 
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where 



lows 



jf =71 -2 l2 e°X + 73 e z °X 
4 V) =71-72 (e- 2a + e°)X + 73 e~° X 2 , (6) 
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and 



71 = 3 + 3 Qf3 + Q!4 , 72 = 1 + 2 0:3 + Q!4 , 73 = «3 + Qf4 

H = a/(aN) , H f = a/ {an) , S = &/N , 

X = a/a , r = an/(aN). (7) 

We note that, in the isotropic limit (er, S — > 0), we have 
jj^' = J^' , so that the Stiickelberg equation of motion, 
Eq. ([5]), at leading order, gives 



7i- 2 72 X + 7 3 X 2 ~0, 



(8) 



that is X — > constant, which corresponds to the FLRW 
result found in [5(. In the same limit, we can also see 
that H — > constant, as expected. 

Even modes. Let us now consider the perturbations 
which transform as 2D scalars under a spatial rotation 
in the y—z plane (also referred as even modes). Then, 
the perturbed metric for the even sector can be written 
as 

ds 2 = -N 2 (l + 2<$>)dt 2 + 2aNdt[e 2a d x xdx + e^diBdy*] 
+ a 2 e 4a {l + 4>)dx 2 + 2a 2 e' J d x d l (3dxdy l 
+ a 2 e- 2CT [%(l + r) + did, E]dy l dy> , (9) 

while the even-type perturbations of Stiickelberg fields 
read 

(p° = t + tt° , cf) 1 = x + dzir 1 , 4> l = y l + <9V . (10) 
We can then define gauge invariant combinations as fol- 
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The first four definitions do not refer to the Stiickelberg 



perturbations and are thus already present in GR [11 
However, the additional three degrees arise from the 
breaking of general coordinate invariance by the non zero 
expectation value of the Stiickelberg fields. 

In order to find the behavior of the perturbations, we 
proceed as usual by expanding the action at second order 
in the perturbation fields, then by employing the Fourier 
plane-wave decompositions, as in cxp[i(fcia; + kiy 1 )]. The 
degrees of freedom arising from the perturbations, 
namely 4>, B and x, are nondynamical, thus can be in- 
tegrated out. Furthermore, the kinetic term for the f„ 
is proportional to the background equations of motion, 
so that this degree of freedom is also nondynamical. We 
interpret this field as the would-be BD ghost, which is 
eliminated in this theory by construction. 

In the massless theory (i.e. GR), using the constraint 
equations also removes the degrees j3 n , E n , leaving only 
i\) in the action, which becomes one of two gravity wave 
polarizations in the isotropic limit. However, in our case, 
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due to the nonzero mass of the graviton, these two de- 
grees of freedom are dynamical, in general. 

Thus, the Lagrangian for even- type perturbations in 
vacuum has three physical propagating modes, V a , (a = 
1,2,3). Assuming small deviation from FLRW, with 
\a\ <C 1 and \S/H\ < 1, we study the kinetic matrix 
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Thanks to the 2D rotational symmetry on the y—z plane, 
the action depends on kx = \Jk\ + fc 2 , instead of the in- 
dividual components. The eigenvalues of /C a f>, at leading 
order in small anisotropy expansion, are 
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where we defined Mg W 
introduced the physical momenta 



m2(l-r)X 2 ( 72 - 73 X), and 
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scales. Therefore, at sufficiently large scales, this mode 
cannot be integrated out from the low energy effective 
theory. This feature makes the FLRW background unsta- 
ble for massive gravity. As a consequence, the homoge- 
neous and isotropic cosmology cannot be accommodated 
in the nonlinear massive gravity theory. 

Odd modes. Let us now discuss the odd sector (i.e. 
the divergence-less part of the modes which transform 
as 2D vectors under a rotation in the y—z plane). The 
perturbed metric we consider is 

ds 2 = -N 2 dt 2 + 2ae~ cr Nv l dtdy z + 2a 2 e' J d x X i dxdy i 

+ a 2 e ia dx 2 + a 2 e- 2,J {5 l] + ,1, , /.dy , (17) 

where duhj) = (dihj +dihj)/2 and d l Vi = d l \ = d l hi = 
0. For the Stuckelberg fields, we consider instead 



t , (f) 1 = x ,4? = y l + tt' 



(18) 



where dirf = 0. Since the vectors are defined on the 2D 
y—z plane, the transverse condition can be used to reduce 
each of these vectors to a single degree of freedom 



The kinetic term k± which is the only eigenvalue that does 
not vanish in isotropic limit, corresponds to one of the 
gravity wave polarizations in FLRW. Once small but non- 
vanishing anisotropy is introduced, two additional even 
modes acquire nonzero kinetic terms at quadratic order. 
More importantly, from t|l 3|) , we see that K2 and K3 have 
opposite signs, regardless of the parameters of the the- 
ory. Thus, we conclude that in the isotropic limit, one 
of the new degrees is always a ghost. Assuming that 
a(l — r) > (which turns out to be the condition for 
stability in the odd sector, as we show later), the ghost 
mode is associated with the eigenvalue K2 < 0. 

We conclude the discussion of the even modes by pre- 
senting their dispersion relations. We first make a field 
redefinition into new field basis fields W a defined such 
that the kinetic action can be written as 
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Ndtdk L d 2 k T a 3 (^Vab^- J , (15) 



where rj a b = diag(l, — 1, 1). The mass spectrum can 
be determined either by studying the equation for the 
frequency-discriminant, or equivalcntly, by performing a 
Lorentz transformation to diagonalize the frequency ma- 
trix. Eventually, we find 
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with 



; 2 w| < in general, and lu 2 < by assuming 
cr(l — r) > 0. We note that the dispersion relation cor- 
responding to the ghost, u>2, becomes smaller at larger 



Vi = e^djV, Aj = e^djX, h t = e/djh, ^ = e/ djir odd , 

where is a unit anti-symmetric tensor with e 2 3 = 
—e 2 = 1. Also for the odd modes we can introduce 
gauge invariant combinations as follows 

-la 



a e 



2N 

,-3<T 



h, 



A = A - - 

tin = 7T dd 



-h. 

2 



(19) 



Using these fields, the second-order resulting action de- 
pends on the three perturbations (v, A, h w ). Among 
these, v does not have any time derivatives and can be 
removed by solving its own constraint equation. In Gen- 
eral Relativity, this operation also removes h n and the 
final action can be written in terms of A only. However, 
in this nonlinear theory of massive gravity, we expect the 
field hn to remain in the action as an extra degree of 
freedom coming from the Stuckelberg sector. 
After a further field redefinition, 
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the quadratic action, for small anisotropy, takes the fol- 
lowing form 
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where the two modes decouple at leading order in the 
small anisotropy expansion, with coefficients 
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and ct 



<i<i 



(1 — r 2 )/(2er). Thus, at leading order, we 



identify the mode Qi with one of the gravity wave polar- 
izations in the FLRW background The extra degree 
of freedom Q 2 is massless and has sound speed c 0( jd- In 
order for this mode to be stable, we require the kinetic 
term for Q 2 to be positive, that is 



(1 - r) a > . 



(23) 



In this case, also 



xld 



becomes positive, and the odd 



mode Q 2 is, in general, free from ghost instabilities. 

Conclusions. In the search for a theory which could 
explain the dark energy enigma, the nonlinear massive 
gravity, recently introduced in |l|, nas raised lots of in- 
terest among both theoretical and experimental physi- 
cists, thanks to its implications in our understanding of 
fundamental forces, if it is theoretically consistent and 
observationally viable. 

This theory admits homogeneous and isotropic solu- 
tions, and it has been shown in [H[ that, out of the five 
modes which would be typically expected in this the- 
ory, only two actually propagate at the linearized level. 
Therefore, it is of interest to investigate the reason for 
this unexpected feature. 

We propose here that this phenomenon is due to the 
high symmetry structure of the FLRW background. Ac- 
cordingly, we have studied the small anisotropy limit of 
the Bianchi-I manifold and found that there is always a 
ghost mode in the even sector, with a propagation speed 
that diverges in the isotropic limit. Furthermore, this 
mode does not have a mass gap; its frequency tends to 
zero for small values of the momentum. Therefore, at 
sufficiently large scales, the frequency cannot be consid- 
ered as large compared to the ultraviolet cutoff of the 
theory. As a consequence, the ghost mode cannot be 
integrated out in general, and the almost-isotropic back- 
ground becomes unstable under production of negative 
energy quanta. 

Although our analysis is linear, the terms in the 
quadratic action with coefficients proportional to the 
small anisotropy can be interpreted as the leading order 
nonlinear perturbations in a pure FLRW universe. The 
presence of a mode with negative kinetic term indicates 
that a homogeneous and isotropic universe in nonlinear 
massive gravity is unstable. 

Our conclusion about ghost instability is far more 
general than it appears d, despite the simplicity of the 



analysis presented above. This is because, whenever a 
quadratic kinetic term vanishes, the leading kinetic term 
is generically cubic and thus can easily become negative, 
signaling the existence of ghost at the nonlinear level. 
Moreover, the other type of homogeneous and isotropic 
solutions (in the non-self-accelerating branch) suffer from 



13] 



as 



ghost instability already at the linearized level 
expected from the classical work of Higuchi [l4| • There- 
fore, all homogeneous and isotropic backgrounds, as well 
as most (if not all) of known spherically-symmetric inho- 
mogencous solutions, are unstable. 
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